### Expdata
fit = glm(now ~ larea, family = binomial(), data = expdata)
a=fit$coef[1]
b=fit$coef[2]
plot(expdata$larea, expdata$now)
curve(p.fun,0,10,add=T)

### Coldata

fitC = glm(now ~ larea, family = binomial(), data = coldata)
a = fitC$coef[1]
b = fitC$coef[2]
plot(coldata$larea, coldata$now)
curve(p.fun,0,10,add=T)

summary(fitC)